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ticular, that the complete account of scattering processes leads to a renormalization 
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with a thermodynamic Lyapunov function. 
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Introduction 



In this paper, explicit method of constructing approximations (the Triangle 
Entropy Method [1,2]) is developed for nonequilibrium problems of Boltz- 
mann's-type kinetics, i.e. when the standard moment variables become insuf- 
ficient. This method enables one to treat any complicated nonlinear functionals 
that fit best the physics of a problem (such as, for example, rates of processes) 
as new independent variables. 

The method is applied to the problem of derivation of hydrodynamics from the 
Boltzmann equation. New macroscopic variables are introduced (moments of 
the Boltzmann collision integral, or scattering rates). They are treated as in- 
dependent variables rather than as infinite moment series. This approach gives 
the complete account of rates of scattering processes. Transport equations for 
scattering rates are obtained (the second hydrodynamic chain), similar to the 
usual moment chain (the first hydrodynamic chain). 

Using the triangle entropy method, three different types of the macroscopic 
description are considered. The first type involves only moments of distribution 
functions, and results coincide with those of the Grad method in the Maximum 
Entropy version. The second type of description involves only scattering rates. 
Finally, the third type involves both the moments and the scattering rates (the 
mixed description). 

The second and the mixed hydrodynamics are sensitive to the choice of the 
collision model. The second hydrodynamics is equivalent to the first hydrody- 
namics only for Maxwell molecules, and the mixed hydrodynamics exists for 
all types of collision models excluding Maxwell molecules. 

Various examples of the closure of the first, of the second, and of the mixed 
hydrodynamic chains are considered for the hard spheres model. It is shown, 
in particular, that the complete account of scattering processes leads to a 
renormalization of transport coefficients. 

The method gives the explicit solution for the closure problem, provides ther- 
modynamic properties of reduced models, and can be applied to any ki- 
netic equation with a thermodynamic Lyapunov function, for example, to the 
Fokker-Planck Equation. 

Reduction of description for dissipative kinetics assumes (explicitly or implic- 
itly) the following picture (Fig. la): There exists a manifold of slow motions 
Q s i ow in the space of distributions. From the initial conditions the system goes 
quickly in a small neighborhood of the manifold, and after that moves slowly 
along it. The manifold of slow motion (slow manifold, for short) must be pos- 
itively invariant: if a motion starts on the manifold at t , then it stays on the 
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manifold at t > to. In some neighborhood of the slow manifold the directions of 
fast motion could be defined. Of course, we always deal not with the invariant 
slow manifold, but with some approximate (ansatz) slow manifold Q. 

There are three basic problems in the model reduction: 

(1) How to construct an approximate slow manifold; 

(2) How to project the initial equation onto the constructed approximate 
slow manifold, i.e. how to split motions into fast and slow; 

(3) How to improve the constructed manifold and the projector in order to 
make the manifold more invariant and the motion along it slower. 

The first problem is often named "the closure problem" , and its solution is the 
closure assumption; the second problem is "the projection problem". Some- 
times these problems are discussed and solved simultaneously (for example, for 
the quasiequilibrium, or, which is the same, for MaxEnt closure assumptions 
[3,4,5,6,7,9]). Sometimes solution of the projection problem after construction 
of ansatz takes a long time. The known case of such a problem gives us the 
Tamm-Mott-Smith approximation in the theory of shock waves (see, for ex- 
ample, [10]). However if one has constructed the closure assumption which is 
at the same time the invariant manifold [10,12,13], then the projection prob- 
lem disappears, because the vector field is always tangent to the invariant 
manifold. In this paper, we would like to add several new tools to the col- 
lection of methods for solution the closure problem. The second problem was 
discussed in Ref. [11]. We do not discuss here the third main problem of model 
reduction: How to improve the constructed manifold and the projector in or- 
der to make the manifold more invariant and the motion along it more slow. 
This discussion can be found in various works [10,12,13,14,15], and the broad 
review of the methods for invariant manifolds construction was performed in 
Refs. [16,17]. 

Our standard example in this paper is the Boltzmann equation, but most of 
the methods can be applied to an almost arbitrary kinetic equation with a 
convex thermodynamic lyapunov function. Let us discuss the initial kinetic 
equation as an abstract ordinary differential equation 1 , 

1 Many of partial differential kinetic equations or integro-differential kinetic equa- 
tions with suitable boundary conditions (or conditions at infinity) can be discussed 
as abstract ordinary differential equation in appropriate space of functions. The 
corresponding semigroup of shifts in time can be considered too. Sometimes, when 
an essential theorem of existence and uniqueness of solution is not proven, it is 
possible to discuss a corresponding shift in time with the support of physical sense: 
the shift in time for physical system should exist. Benefits from the latter approach 
are obvious as well as its risk. 
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Fig. 1. a) Fast-slow decomposition. Bold dashed line - slow invariant manifold; 
bold line - approximate invariant manifold; several trajectories and correspondent 
directions of fast motion are presented schematically, b) The geometrical structures 
of model reduction: U is the phase space, J(f) is the vector field of the system under 
consideration: df/dt = J{f), f2 is an ansatz manifold, Tf is the tangent space to the 
manifold Q at the point /, PJ{f) is the projection of the vector J(f) onto tangent 
space Tf, A = (1 — P)J(f) is the defect of invariance, the affine subspace / + ker P 
is the plain of fast motions, and A G ker P. 

where / = f(q) is the distribution function, q is the point in particle phase 
space (for the Boltzmann equation), or in configuration space (for the Fokker- 
Planck equation). This equation is denned in some domain U of a vector space 
of admissible distributions E. 

The dissipation properties of the system (1) are described by specifying the 
entropy S, the distinguished Lyapunov function which monotonically increases 
along solutions of equation (1). In a certain sense, this Lyapunov function is 
more fundamental than the system (1) itself. That is, usually, the entropy is 
known much better than the right hand side of equation (1). 

We assume that a concave functional S is defined in U, such that it takes 
maximum in an inner point /* G U. This point is termed the equilibrium. 

For any dissipative system (1) under consideration in U, the derivative of S 
due to equation (1) must be nonnegative, 



where DfS is the linear functional, the differential of the entropy. 

For dissipative systems, we always keep in mind the following picture (Fig. lb). 
The vector field J(f) generates the motion on the phase space U: df/dt = 
</(/). An ansatz manifold Q is given, it is the current approximation to the 
invariant manifold. 

The projected vector field PJ(f) belongs to the tangent space Tf, and the 
equation df/dt = PJ(f) describes the motion along the ansatz manifold Q 
(if the initial state belongs to Q). 



dS_ 



(D f S)(J(f))>0, 



(2) 
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Fig. 2. Parametrization by macroscopic variables: linear (a), nonlinear (b) and 
layer-linear (c). Thin arrows illustrate the bijection M <-> /m- a) Moment 
parametrization in fast-slow decomposition. Dashed lines - the plains of constant 
value of moments M. These plains coincide with directions of fast motion in the 
moment approximation, b) Non-linear macroscopic parametrization in fast-slow de- 
composition. Dashed curves - the surfaces of constant value of macroscopic variables 
M. Plains of fast motion are tangent to these surfaces, c) Nonlinear, but layer-linear 
macroscopic parametrization in fast-slow decomposition. The surfaces of constant 
value of macroscopic variables M (dashed lines) are plain, but the dependence m(f) 
is nonlinear. Plains of fast motion coincide with these plains. 



The choice of the projector P might be very important. There is some "dual- 
ity" between accuracy of slow invariant manifold approximation and restric- 
tions on the projector choice. If Q is an exactly invariant manifold, then the 
vector field J(f) is tangent to f2, and all projectors give the same result. If 
Q gives a good smooth approximation for such an invariant manifold, then 
the set of admissible projectors is rather broad. On the other hand, there is 
the unique choice of the projector applicable for every (arbitrary) ansatz Q 
[10,11], any other choice leads to violation of the Second Law in projected 
equations. 

In the initial geometry of the fast-slow decomposition (Figs, la, lb) the "slow 
variables" (or "macroscopic variables") are internal coordinates on the slow 
manifold, or on its approximation Q. It is impossible, in general, to define these 
macroscopic variables as functionals of / outside these manifolds. Moreover, 
this definition cannot be unique. 

The moment parametrization starts not from the manifold, but from the 
macroscopic variables defined in the whole U (Fig. 2a), and for the given vari- 
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ables it is necessary to find the corresponding slow manifold. Usually, these 
slow variables are linear functions (functionals), for example, hydrodynamic 
fields (density, momentum density, and pressure) are moments of one-particle 
distribution functions f(x,v). The moment vector M is the value of linear 
operator m: M = m(f). The moments value serve as internal coordinates on 
the (hypothetic) approximate slow manifold Q. It means that points of Q are 
parameterized by M, fl = {/m}, and the consistency condition holds: 

m(f M ) = M. 



In the example with one-particle function / and hydrodynamic fields m(f) 
it means that slow manifold consists of distribution f(x,v) parameterized by 
their hydrodynamic fields (as in the Chapman-Enskog method [34]). For a 
given fl = {/m}, the moment equation has a very simple form: 

^-=m(J(f M )), (3) 



and the correspondent equation for the projected motion on the manifold 
^ = {I'm} is 

^ = (D M f M )m(J(f M )), (4) 



Where D M f M is differential of the parametrization M h- > fu- 

How to find a manifold Q = {Jm} for a given moment parametrization? A good 
initial approximation is the quasiequilibrium (or MaxEnt) approximations. 
The basic idea is: in the fast motion the entropy should increase, hence, the 
point of entropy maximum on the plane of rapid motion is not far from the 
slow manifold (Fig. la). If our moments M are really slow variables, and don't 
change significantly in rapid motion, then the manifold of conditional entropy 
maxima fu'- 

S(f) - max, m(f) = M (5) 



could serve as appropriate ansatz for slow manifold. 



Most of the works on nonequilibrium thermodynamics deal with quasiequilib- 
rium approximations and corrections to them, or with applications of these 
approximations (with or without corrections). This viewpoint is not the only 
possible but it proves very efficient for the construction of a variety of useful 
models, approximations and equations, as well as methods to solve them. From 
time to time it is discussed in the literature, who was the first to introduce 
the quasiequilibrium approximations, and how to interpret them. At least a 
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part of the discussion is due to a different role the quasiequilibrium plays in 
the entropy- conserving and the dissipative dynamics. The very first use of the 
entropy maximization dates back to the classical work of G. W. Gibbs [19], 
but it was first claimed for a principle of informational statistical thermo- 
dynamics by E. T. Jaynes [3]. Probably the first explicit and systematic use 
of quasiequilibria to derive dissipation from entropy-conserving systems was 
undertaken by D. N. Zubarev. Recent detailed exposition is given in [4]. For 
dissipative systems, the use of the quasiequilibrium to reduce description can 
be traced to the works of H. Grad on the Boltzmann equation [20] . A review 
of the informational statistical thermodynamics was presented in [21]. The 
connection between entropy maximization and (nonlinear) Onsager relations 
was also studied [22,23]. The viewpoint of the present authors was influenced 
by the papers by L. I. Rozonoer and co-workers, in particular, [5,6,7]. A de- 
tailed exposition of the quasiequilibrium approximation for Markov chains is 
given in the book [18] (Chap. 3, Quasiequilibrium and entropy maximum, pp. 
92-122), and for the BBGKY hierarchy in the paper [9]. The maximum en- 
tropy principle was applied to the description the universal dependence the 
three-particle distribution function F 3 on the two-particle distribution func- 
tion F 2 in classical systems with binary interactions [24]. For a discussion the 
quasiequilibrium moment closure hierarchies for the Boltzmann equation [6] 
see the papers [1,2,25]. A very general discussion of the maximum entropy 
principle with applications to dissipative kinetics is given in the review [26]. 
Recently the quasiequilibrium approximation with some further correction was 
applied to description of rheology of polymer solutions [27,28] and of ferrofluids 
[29,30]. Quasiequilibrium approximations for quantum systems in the Wigner 
representation [31,32] was discussed very recently [33]. 

Formally, for quasiequilibrium approximation the linearity of the map / i— > 
m(f) is not necessary, and the optimization problem (5) could be studied for 
nonlinear conditions m(f) = M (Fig. 2b). Nevertheless, the problem (5) with 
nonlinear conditions loose many important properties caused by concavity of 
S. The technical compromise is the problem with nonlinear map to, but linear 
restrictions m(f) = M. It is possible when preimages of points for the map to 
are plain (Fig. 2c). Such a "layer-linear" approximation for a generic smooth 
map to : / i— > M could be created as follows. Let Qq be a smooth submanifold 
in U. In some vicinity of Q we define a map to 

m(f) = to (/o) if (Dm ) fo (f - fo) = 0, (6) 

where fo are points from Q and (_Dmo)/ is the differential of m at the point 
fo- This definition means that m(f) = m (fo) if m Q (f) coincides with too(/o) 
in linear approximation. 

The layer-linear parametrization was introduced in Ref. [35] for construction 
of generalized model equations for the Boltzmann kinetics. 
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Eq. (6) defines a smooth layer-linear map m in a vicinity of VLq under some 
general transversality condition. 

Let an initial approximation VLq for the slow manifold be unsatisfactory. Two 
basic ways for improvement are: (i) manifold correction and (ii) manifold ex- 
tension. On the first way we should find a shifted manifold that is better ap- 
proximate slow invariant manifold. The list of macroscopic variables remains 
the same. On the second way we extend the list of macroscopic variables, and, 
hence, extend the manifold Vt. The Chapman-Enskog method [34] gives the 
example of manifold correction in the form of Taylor series expansion, the 
direct Newton method could give better results [10,12,16,17,39]. The second 
way is the essence of EIT - extended irreversible thermodynamics [44]. This 
paper is focused on the manifold extensions also. 

Usually moments are graduated in a natural order, by degree of polynomials: 
concentration (zero order of velocity), average momentum density (first order), 
kinetic energy (second order), stress tensor (second order), heat flux (third or- 
der), etc. Normal logic of EIT is extension of the list of variables by addition 
of the next irreducible moment tensor. But there is another logic. In general, 
for the set of moments M that parametrizes Qq a time derivative is a known 
function of /: dM/dt = Fm(/)- We propose to construct new moments from 
F M (f). It allows to achieve a best possible approximation for dM/dt through 
extended variables. For this nonlinear variables we use the layer-linear ap- 
proximation (6), ass well, as a layer-quadratic approximation for the entropy. 
This (layer) linearization of the problem near current approximation follows 
lessons of the Newton method. 

It should be stressed that "layer-linear" does not mean "linear" , and the mod- 
ified choice of new variables implies no additional restrictions, but it pretends 
to be just a more direct way to dynamic invariance. Below this approach is 
demonstrated on the Boltzmann equation. 



1 Difficulties of classical methods of the Boltzmann equation the- 
ory 

The Boltzmann equation remains the most inspiring source for the model 
reduction problems. The first systematic and (at least partially) successful 
method of constructing invariant manifolds for dissipative systems was the 
celebrated Chapman-Enskog method [34] for the Boltzmann kinetic equation. 
The main difficulty of the Chapman-Enskog method [34] are "nonphysical" 
properties of high-order approximations. This was stated by a number of au- 
thors and was discussed in detail in [36]. In particular, as it was noted in [37], 
the Burnett approximation results in a short-wave instability of the acoustic 
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spectra. This fact contradicts the fi-theorem (cf. in [37]). The Hilbert expan- 
sion contains secular terms [36]. The latter contradicts the if -theorem. 

The other difficulties of both of these methods are: the restriction upon the 
choice of the initial approximation (the local equilibrium approximation), the 
requirement for a small parameter, and the usage of slowly converging Taylor 
expansion. These difficulties never allow a direct transfer of these methods on 
essentially nonequilibrium situations. 

The main difficulty of the Grad method [20] is the uncontrollability of the 
chosen approximation. An extension of the list of moments can result in a 
certain success, but it can also give nothing. Difficulties of moment expansion 
in the problems of shock waves and sound propagation are discussed in [36]. 

Many attempts were made to refine these methods. For the Chapman-Enskog 
and Hilbert methods these attempts are based in general on some better re- 
arrangement of expansions (e.g. neglecting high-order derivatives [36], reex- 
panding [36], Pade approximations and partial summing [1,38,40,41], etc.). 
This type of work with formal series is wide spread in physics. Sometimes 
the results are surprisingly good - from the renormalization theory in quan- 
tum fields to the Percus-Yevick equation and the ring-operator in statistical 
mechanics. However, one should realize that success cannot be guaranteed. 
Moreover, rearrangements never remove the restriction upon the choice of the 
initial local equilibrium approximation. 

Attempts to improve the Grad method are based on quasiequilibrium approx- 
imations [5,6]. It was found in [6] that the Grad distributions are linearized 
versions of appropriate quasiequilibrium approximations (see also [1,2,25]). A 
method which treats fluxes (e.g. moments with respect to collision integrals) 
as independent variables in a quasiequilibrium description was introduced in 
[1,2,42,43], and will be discussed later. 

The important feature of quasiequilibrium approximations is that they are 
always thermodynamic, i.e. they are consistent with the if -theorem by con- 
struction. 



2 Triangle Entropy Method 

In the present subsection, which is of introductory character, we shall refer, to 
be specific, to the Boltzmann kinetic equation for a one-component gas whose 
state (in the microscopic sense) is described by the one-particle distribution 
function f(v,x,t) depending on the velocity vector v = {vk} 3 k=1 , the spatial 
position x = {x k } k=1 and time t. The Boltzmann equation describes the 
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evolution of / and in the absence of external forces is 
d t f + v k d k f = Q(f,f), 



(7) 



where d t = d/dt is the time partial derivative, = d/dx^ is partial derivative 
with respect to fc-th component of x, summation in two repeating indices is 
assumed, and Q(f,f) is the collision integral (its concrete form is of no im- 
portance right now, just note that it is functional-integral operator quadratic 
with respect to /). 

The Boltzmann equation possesses two properties principal for the subsequent 
reasoning. 

1. There exist five functions if) a (v) (additive collision invariants), 

l,v,v 2 

such that for any their linear combination with coefficients depending on x, t 
and for arbitrary / the following equality is true: 

/ j2a«MMv)Q(f,f)dv = 0, (8) 

J a=l 

provided the integrals exist. 

2. The equation (7) possesses global Lyapunov functional: the if-function, 

H{t) = H[f] = J f(v, x, t) lnf(v, x, t) dv dx, (9) 



the derivative of which by virtue of the equation (7) is non-positive under 
appropriate boundary conditions: 

dH(t)/dt < 0. (10) 



Grad's method [20] and its variants construct closed systems of equations 
for macroscopic variables when the latter are represented by moments (or, 
more general, linear functionals) of the distribution function / (hence their 
alternative name is the "moment methods"). The entropy maximum method 
for the Boltzmann equation is of particular importance for the subsequent 
reasoning. It consists in the following. A finite set of moments describing the 
macroscopic state is chosen. The distribution function of the quasiequilibrium 
state under given values of the chosen moments is determined, i.e. the problem 
is solved 

H[f] - min, for M,[/] = M h i — l,...,k, (11) 
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where Mi[f] are linear functionals with respect to /; Mj are the corresponding 
values of chosen set of k macroscopic variables. The quasiequilibrium distri- 
bution function f*(v, M(x, t)), M = {Mi, . . . , M&}, parametrically depends 
on Mj, its dependence on space x and on time t being represented only by 
M(x, t). Then the obtained /* is substituted into the Boltzmann equation (7), 
and operators M are applied on the latter formal expression. 

In the result we have closed systems of equations with respect to Mi(x,t), 
i = 1, . . . , k: 

dtM, + M t [v k d k f*(v, M)] = Mi[Q(f*(v, M), />, M))]. (12) 

The following heuristic explanation can be given to the entropy method. A 
state of the gas can be described by a finite set of moments on some time scale 
9 only if all the other moments ("fast") relax on a shorter time scale time 
r, r << 9, to their values determined by the chosen set of "slow" moments, 
while the slow ones almost do not change appreciably on the time scale r. 
In the process of the fast relaxation the //-function decreases, and in the 
end of this fast relaxation process a quasiequilibrium state sets in with the 
distribution function being the solution of the problem (11). Then "slow" 
moments relax to the equilibrium state by virtue of (12). 

The entropy method has a number of advantages in comparison with the classi- 
cal Grad's method. First, being not necessarily restricted to any specific system 
of orthogonal polynomials, and leading to solving an optimization problem, it 
is more convenient from the technical point of view. Second, and ever more 
important, the resulting quasiequilibrium //-function, H*(M) = H[f*(v, M)], 
decreases due of the moment equations (12). 

Let us note one common disadvantage of all the moment methods, and, in 
particular, of the entropy method. Macroscopic parameters, for which these 
methods enable to obtain closed systems, must be moments of the distribution 
function. On the other hand, it is easy to find examples when the interesting 
macroscopic parameters are nonlinear functionals of the distribution function. 
In the case of the one-component gas these are the integrals of velocity poly- 
nomials with respect to the collision integral Q(f, f) of (7) (scattering rates 
of moments). For chemically reacting mixtures these are the reaction rates, 
and so on. If the characteristic relaxation time of such nonlinear macroscopic 
parameters is comparable with that of the "slow" moments, then they should 
be also included into the list of "slow" variables on the same footing. 

In this paper for constructing closed systems of equations for non-linear (in a 
general case) macroscopic variables the triangle entropy method is used. Let 
us outline the scheme of this method. 
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Let a set of macroscopic variables be chosen: linear functional M[f] and 
nonlinear functionals (in a general case) N[f]: 

M[f] = {M,[fl ...,M k [f]}, N[f] = {#![/], . . . , N[f]} . 

Then, just as for the problem (11), the first quasiequilibrium approximation 
is constructed under fixed values of the linear macroscopic parameters M: 

H[f] -> min for M t [f\ = M t , i = 1, . . . , k, (13) 



and the resulting distribution function is f*(v,M). After that, we seek the 
true quasiequilibrium distribution function in the form, 

/ = /•(! + ¥>), (14) 



where <p is a deviation from the first quasiequilibrium approximation. In order 
to determine tp, the second quasiequilibrium approximation is constructed. Let 
us denote AH[f*, ip] as the quadratic term in the expansion of the if -function 
into powers of ip in the neighbourhood of the first quasiequilibrium state /*. 
The distribution function of the second quasiequilibrium approximation is the 
solution to the problem, 



AH[f*,(p] -> min for 
Mi[/>] = 0, i = l,...,k, 

AN j \r,<p]=AN j , 3 = 1,..., I, (15) 



where ANj are linear operators characterizing the linear with respect to ip 
deviation of (nonlinear) macroscopic parameters Nj from their values, N* = 
Nj[f*], in the first quasiequilibrium state. Note the importance of the homo- 
geneous constraints Mi[f*(p] = in the problem (15). Physically, it means 
that the variables ANj are "slow" in the same sense, as the variables Mj, at 
least in the small neighborhood of the first quasiequilibrium /*. The obtained 
distribution function, 

f = f*(v,M)(l + <p**(v,M,AN)) (16) 



is used to construct the closed system of equations for the macroparameters 
M, and AN. Because the functional in the problem (15) is quadratic, and all 
constraints in this problem are linear, it is always explicitly solvable. 

Further in this section some examples of using the triangle entropy method 
for the one-component gas are considered. Applications to chemically reacting 
mixtures were discussed in [42]. 
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3 Linear macroscopic variables 



Let us consider the simplest example of using the triangle entropy method, 
when all the macroscopic variables of the first and of the second quasiequilib- 
rium states are the moments of the distribution function. 



3. 1 Quasiequilibrium projector 

Let p,\(v), . . . , Hk( v ) be the microscopic densities of the moments 

M 1 {x,t),...,M k {x,t) 
which determine the first quasiequilibrium state, 

Mi(x, t)= J m(v)f(v, x, t) dv, (17) 
and let vi(v), . . . , vi{v) be the microscopic densities of the moments 

determining together with (7) the second quasiequilibrium state, 

Ni(x, t) = J Vi(v)f(v, x, t) dv. (18) 



The choice of the set of the moments of the first and second quasiequilibrium 
approximations depends on a specific problem. Further on we assume that 
the microscopic density /i = 1 corresponding to the normalization condition is 
always included in the list of microscopic densities of the moments of the first 
quasiequilibrium state. The distribution function of the first quasiequilibrium 
state results from solving the optimization problem, 

H[f\ = J f(v)\nf(v)dv^mm (19) 



for 



J pii(v)f(v) dv — Mi,i = !,...,&;. 



Let us denote by M = {Mi, . . . , M^} the moments of the first quasiequilibrium 
state, and by f*(v, M) let us denote the solution of the problem (19). 
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The distribution function of the second quasiequilibrium state is sought in the 
form, 

f = r(v,M)(l + V>)- (20) 



Expanding the if- function (9) in the neighbourhood of f*(v, M) into powers 
of ip to second order we obtain, 

AH(x,t) = AH[f*,ip]=H*(M) + J f*(v,M) \nf*(v, M)<p(v) dv 

+ 1 -Jr(v,M) V 2 (v)dv, (21) 

where H*(M) = H[f*(v,M)] is the value of the iT-function in the first 
quasiequilibrium state. 

When searching for the second quasiequilibrium state, it is necessary that the 
true values of the moments M coincide with their values in the first quasiequi- 
librium state, i.e., 

Mi = J ^(«)/>,M)(l + dt; 

= J K(v)f*(v, M) dv = M*, i = 1, . . . , k. (22) 

In other words, the set of the homogeneous conditions on (p in the problem 
(15), 

j m{v)r{v, M) V {v) dv = 0, i = 1, . . . , k, (23) 

ensures a shift (change) of the first quasiequilibrium state only due to the new 
moments JVi, . . ., JVj. In order to take this condition into account automati- 
cally, let us introduce the following structure of a Hilbert space: 

(1) Define the scalar product 

tyiM = J f*(v,M)Mv)Mv)dv. (24) 

(2) Let be the linear hull of the set of moment densities 

{/i! («),..., M«)}. 

Let us construct a basis of {ei(v), . . . , e r (v)} that is orthonormal in 
the sense of the scalar product (24): 

(ei,ej) = 5ij, (25) 
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i,j — l,...,r; Sij is the Kronecker delta. 
(3) Define a projector P* on the first quasiequilibrium state, 

P> = 5>( ei ,# (26) 
i=i 

The projector P* is orthogonal: for any pair of functions ipi,ip2, 

(P> 1 ,(i-P*)^ 2 ) = 0, (27) 

where 1 is the unit operator. Then the condition (23) amounts to 

P*<p = 0, (28) 

and the expression for the quadratic part of the if-function (21) takes 
the form, 

AH[rM = H*(M) + (In P» + (l/2)(^). (29) 

Now, let us note that the function In /* is invariant with respect to the action 
of the projector P*: 

P*m/* = m/*. (30) 

This follows directly from the solution of the problem (19) using of the method 
of Lagrange multipliers: 



/* = exp£A i (M) Aii (tO, 

i=i 

where Aj(M) are Lagrange multipliers. Thus, if the condition (28) is satisfied, 
then from (27) and (30) it follows that 



(lnr,v?) = (P*lnr,(i-P» = 0. 
Condition (28) is satisfied automatically, if AiVj are taken as follows: 

AN i = ((l-P*)v i ,<p),i = l,...,l. (31) 

Thus, the problem (15) of finding the second quasiequilibrium state reduces 
to 

Aif [/«>] - H*(M) = (1/2)^,^) min for 

((i-F*)v i ,<p)=AN i , i = l,...,l. (32) 
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Note that it is not ultimatively necessary to introduce the structure of the 
Hilbert space. Moreover that may be impossible, since the "distribution func- 
tion" and the "microscopic moment densities" are, strictly speaking, elements 
of different (conjugate one to another) spaces, which may be not reflexive. 
However, in the examples considered below the mentioned difference is not 
manifested. 

In the remainder of this section we demonstrate how the triangle entropy 
method is related to Grad's moment method. 

3.2 Ten-moment Grad approximation. 

Let us take the five additive collision invariants as moment densities of the 
first quasiequilibrium state: 

2 

fflV 

Ho = 1; fi k = v k (k = 1,2,3); /x 4 = (33) 

where v k are Cartesian components of the velocity, and m is particle's mass. 
Then the solution to the problem (19) is the local Maxwell distribution func- 
tion /(°)(«,aj,*): 

where 

n(x,t) = J f(v) dv is local number density, 

u(x,t) = n~ 1 (x : t) J f(v)vdv is the local flow density, 

T(x,t) = -^n~ 1 (x,t) J f(v)(v — u(x,t)) 2 dv is the local temperature, 

ks is the Boltzmann constant. 

Orthonormalization of the set of moment densities (33) with the weight (34) 
gives one of the possible orthonormal basis 



5/cb7" — m(v — it) 2 
6o= (10n)V2jfc B T ' 

6k ~ (nk B Ty/i W 
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e 4 



m(v — u) 2 



{lhn) l ' 2 k B T' 



For the moment densities of the second quasiequilibrium state let us take, 
u ik = mviVk, i,k = 1,2,3. (36) 



Then 



(1 - P {0) )v ik = m(vi - Ui)(v k - u k ) - ^S ik m(v - u) 2 , (37) 



and, since ((1 - P (0) )z/ ife , (i - P (0) )z/ is ) = (5u5 ks + 5 k i5 is )Pk B T/m, where P = 
nk B T is the pressure, and a ik = (f, (1 — P^)v ik ) is the traceless part of 
the stress tensor, then from (20), (33), (34), (37) we obtain the distribution 
function of the second quasiequilibrium state in the form 

/ = / (0) (l + [(Vi - Ui)(v h - u k ) - l -5 lk (v - u) 2 ^ (38) 



This is precisely the distribution function of the ten-moment Grad approxi- 
mation (let us recall that here summation in two repeated indices is assumed). 



3.3 Thirteen-moment Grad approximation 



In addition to (33), (36), let us extend the list of moment densities of the 
second quasiequilibrium state with the functions 

c mViv2 ■ 1 o q (on\ 
& = ,2 = 1,2,3. (39) 



The corresponding orthogonal complements to the projection on the first 
quasiequilibrium state are 

(1 - P (0) )6 = f (t* - «0 ((« - u) 2 - ^ . (40) 

The moments corresponding to the densities (1 — P^)£i are the components 
of the heat flux vector qf. 

ft = (v,(i-p (0) )ei)- (4i) 
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Since 



((i-P (0) )&,(i-P (0) K) = o, 



for any i, k, /, then the constraints 



((i - pW)v lk , <p) = a lk , ((1 - P(°))6, <p) = qi 



in the problem (32) are independent, and Lagrange multipliers corresponding 
to & are 

f(— )'*■ (42) 

on V m / 



Finally, taking into account (33), (38), (40), (42), we find the distribution 
function of the second quasiequilibrium state in the form 



/ = /(«) ! + 



2Pk B T 



Ui){v k - u k ) - ^5 ik {v - uf^j 



Pk B T 



[Vi - Ui) 



u 



5k K T 



- 1 



(43) 



which coincides with the thirteen-moment Grad distribution function [20]. 

Let us remark on the thirteen-moment approximation. From (43) it follows 
that for large enough negative values of (i>j — Uj) the thirteen-moment dis- 
tribution function becomes negative. This peculiarity of the thirteen-moment 
approximation is due to the fact that the moment density £j is odd-order 
polynomial of V{. In order to eliminate this difficulty, one may consider from 
the very beginning that in a finite volume the square of velocity of a particle 
does not exceed a certain value v^ x , which is finite owing to the finiteness 
of the total energy, and is such that when changing to infinite volume 
qi — > 0, f^ ax — > oo and qi{vi — Ui){v — u) 2 remains finite. 

On the other hand, the solution to the optimization problem (11) does not 
exist (is not normalizable) , if the highest-order velocity polynomial is odd, as 
it is for the full 13-moment quasiequilibrium. 

Approximation (38) yields AH (29) as follows: 

AH = H^+n^, (44) 
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while AH corresponding to (43) is, 

A U u(0) , °ik®ik QkQkP /,_>. 

AH = HO +n —— + n ——, (45) 



where p = mn, and H^ [ s the local equilibrium value of the if -function 

5 3 3 / 27r \ 

i/ (0) = -nlnn- -nlnP- -n 1 + In — . (46) 



These expressions coincide with the corresponding expansions of the quasiequi- 
librium if -functions obtained by the entropy method, if microscopic moment 
densities of the first quasiequilibrium approximation are chosen as l,i>i, and 
V{Vj, or as 1,V{, ViVj, and ViV 2 . As it was noted in [6], they differs from the 
fi-functions obtained by the Grad method (without the maximum entropy 
hypothesis), and in contrast to the latter they give proper entropy balance 
equations. 

The transition to the closed system of equations for the moments of the first 
and of the second quasiequilibrium approximations is accomplished by pro- 
ceeding from the chain of the Maxwell moment equations, which is equivalent 
to the Boltzmann equation. Substituting / in the form of /^(l + y?) into 
equation (7), and multiplying by //,(«), and integrating over v, we obtain 



d t (l, P (0 Vi(«)) + dt(<p(v),Hi(v)) + dk(v k <p(v),fii(v)) 

+d k (v k , ta(v)) = M Q \pn,(p]. (47) 



Here 



MqK <p\ = J Q(/(°)(l + <p), + <p))H(v) dv 

is a "moment" (corresponding to the microscopic density) Pi(v) with respect 
to the collision integral (further we term Mq the collision moment or the 
scattering rate). Now, if one uses / given by equations (38), and (43) as a 
closure assumption, then the system (47) gives the ten- and thirteen-moment 
Grad equations, respectively, whereas only linear terms in ip should be kept 
when calculating Mq. 

Let us note some limitations of truncating the moment hierarchy (47) by means 
of the quasiequilibrium distribution functions (38) and (43) (or for any other 
closure which depends on the moments of the distribution functions only). 
When such closure is used, it is assumed implicitly that the scattering rates 
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in the right hand side of (47) "rapidly" relax to their values determined by 
"slow" (quasiequilibrium) moments. Scattering rates are, generally speaking, 
independent variables. This peculiarity of the chain (47), resulting from the 
nonlinear character of the Boltzmann equation, distinct it essentially from 
the other hierarchy equations of statistical mechanics (for example, from the 
BBGKY chain which follows from the linear Liouville equation). Thus, equa- 
tions (47) are not closed twice: into the left hand side of the equation for the 
i-th moment enters the (i + l)-th moment, and the right hand side contains 
additional variables - scattering rates. 

A consequent way of closure of (47) should address both sets of variables 
(moments and scattering rates) as independent variables. The triangle entropy 
method enables to do this. 



4 Transport equations for scattering rates in the neighbourhood 
of local equilibrium. Second and mixed hydrodynamic chains 

In this section we derive equations of motion for the scattering rates. It proves 
convenient to use the following form of the collision integral Q(f, /): 

Q(f,f)(v) = J wiv'^v^) (/(tO/M - f(v)f( Vl )) dv'dv'.dv,, (48) 

where v and V\ are velocities of the two colliding particles before the collision, 
v' and v[ are their velocities after the collision, w is a kernel responsible 
for the post-collision relations v'(v, v±) and v[(v, Vi), momentum and energy 
conservation laws are taken into account in w by means of corresponding 5- 
functions. The kernel w has the following symmetry property with respect to 
its arguments: 

w(v[, v'\v, Vi) = w(v[, v'\vi, v) = w{v' , v\ | V\, v) = w(v, Vi | v' \ u' 1 ).(49) 

Let n{v) be the microscopic density of a moment M. The corresponding scat- 
tering rate Mq [/,//] is defined as follows: 

M Q \f,n] = jQ(f,f)(v)^v)dv. (50) 

First, we should obtain transport equations for scattering rates (50), analogous 
to the moment's transport equations. Let us restrict ourselves to the case when 
/ is represented in the form, 

/ = / {0) (i + y?), (si) 
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where is local Maxwell distribution function (34), and all the quadratic 
with respect to tp terms will be neglected below. It is the linear approximation 
around the local equilibrium. 

Since, by detailed balance, 

/ (0) («)/ (0) (»i)=/ (0) («')/ (0) («i) (52) 

for all such (v, v±), (v', v[) which are related to each other by conservation 
laws, we have, 

M Q [/(°),//] = 0, for any /i. (53) 

Further, by virtue of conservation laws, 

M Q [/,P(°V] = 0, for any/. (54) 

From (52)-(54) it follows, 

M [/<°>(1 + <p)M = M Q [p, (I - P<f»)A (55) 
= - f w(v',v[ I v,v 1 )/ (0) (v)/ (0) (v 1 ){(l-P (0) )/i(«)} dw'dwidwidt;. 

We used notation, 

{^j(v)} = ^(v) + - ip(v') - ^K). (56) 

Also, it proves convenient to introduce the microscopic density of the scatter- 
ing rate, /j,q(v): 

li Q {v) = J w(v',v[ I v,v 1 )f(°\v 1 ){(l-P^) f i(v)} dv'dv' 1 dv 1 . (57) 
Then, 

M Q [(p,fi] = -((p,fji Q ), (58) 

where (•, •) is the L 2 scalar product with the weight (34). This is a natural 

scalar product in the space of functions tp (51) (multipliers), and it is obviously 

related to the entropic scalar product in the space of distribution functions 

at the local equilibrium /(°), which is the L 2 scalar product with the weight 
(/ (o))-i. 
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Now, we obtain transport equations for the scattering rates (58). We write 
down the time derivative of the collision integral due to the Boltzmann equa- 
tion, 

d t Q(f, /)(») = TQ(f, f)(v) + RQ(f, f)(v), (59) 



where 



TQ(f, /)(») = J w(v', v[ | v, Vl ) [f(v)v lk d k f( Vl ) + f( Vl )v k d k f(v) 

- f(v')v' lk d k f(v[) - /(«i)Ufc9 fc /(w')] dt/ dt>i dw! dv; (60) 

RQ(f,f)(v) = JwivWilv^) [g(/,/)K)/W) + Q(/,/)K)/K) 

- W,/)(»i)/(tO-Q(/, /)(«)/(»!)] dt/dwidthdw. (61) 

Using the representation, 

9 fe /(°)(t;) = A fc ( V )/( )(t;); (62) 
A k (v) = d k ln(nT" 3 / 2 ) + ^-(v, - Ui )d k u t + ^f^ d k InT, 

and after some simple transformations using the relation 

{A k (v)} = 0, (63) 

in linear with respect to tp deviation from /(°) (51), we obtain in (59): 

TQ(f,f)(v)=d k J w(v',v[ | v,Vi)/ (0) (^i)/ (0) (^)K^)} dv' 1 dv'dv 1 
+ J w(v',v[ | t>,«i)/ (0) («i)/ (0) («) {v k A k (v)} dv'dv'.dv, 
+ J w(v', v[ | v, v 1 )f (0 \v)f(°\v 1 ) [<p(v)A k ( Vl )(v lk - v k ) 

+^{v 1 )A k (v)(v k - v lk ) + ^(v')A k (v[)(v' k - v' lk ) 

+ tp(v' 1 )A k (v')(v > lk -v' k )} dv' 1 dv > dv 1 ; (64) 

RQ(fJ)(v)= J w(v',v\ | v,v 1 )f^(v)f^(v 1 ){av)} dv' 1 dv'dv 1 ; 

£( v ) = J w (v',v\ | w,«i)/ (0) («i)M«)} dw'idw'dwi; (65) 
$<?(/,/)(») (66) 
= «;(«', «i | ^Vi)/ ( °HtO/ (0) («i)M«)} dt/d^d^. 
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Let us use two identities: 



1. From the conservation laws it follows 

{<p(v)} = {(1 - P^Mv)} . (67) 

2. The symmetry property of the kernel w (49) which follows from (49), (52) 

J w(v',v\ I v,v 1 )f^(v 1 )f^\v)g 1 (v) {g 2 (v)} cV dv[ dv, dv (68) 

= J W(V',V\ I »,»l)/ (0) («l)/ (0) («)52(«) {9l(v)} dv'dv'^dv. 

It is valid for any two functions gi, g 2 ensuring existence of the integrals, and 
also using the first identity. 

Now, multiplying (64)-(67) by the microscopic moment density n(v), perform- 
ing integration over v (and using identities (67), (69)) we obtain the required 
transport equation for the scattering rate in the linear neighborhood of the 
local equilibrium: 

-d t AM Q [(p, /i] = -dt(<p,HQ) 
= (v k A k (v),l*Q«l-P (f>) Mv))) 

+d k (y>(v)v k ,n Q ((i-P {0) Mv))) + J w(v',v[ I v,v 1 )f^(v 1 )f^(v) 
x {(1 - P {0) )ii(v)} A k ( Vl )(v lk - v k )<p(v) dv' dv[ dv x dv 

+ (e(«),^((i-p (0) H«)))- (69) 

The chain of equations (69) for scattering rates is a counterpart of the hydro- 
dynamic moment chain (47). Below we call (69) the second chain, and (47) 
- the first chain. Equations of the second chain are coupled in the same way 
as the first one: the last term in the right part of (69) (£,/zq((1 — P^)/x)) 
depends on the whole totality of moments and scattering rates and may be 
treated as a new variable. Therefore, generally speaking, we have an infinite 
sequence of chains of increasingly higher orders. Only in the case of a special 
choice of the collision model - Maxwell potential U = — K,r~ 4 - this sequence 
degenerates: the second and the higher-order chains are equivalent to the first 
(see below). 

Let us restrict our consideration to the first and second hydrodynamic chains. 
Then a deviation from the local equilibrium state and transition to a closed 
macroscopic description may be performed in three different ways for the 
microscopic moment density n(v). First, one can specify the moment M[/j] and 
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perform a closure of the chain (47) by the triangle method given in previous 
subsections. This leads to Grad's moment method. Second, one can specify 
scattering rate Mq \p] and perform a closure of the second hydrodynamic chain 
(69). Finally, one can consider simultaneously both M[/i] and Mq[/i] (mixed 
chain). Quasiequilibrium distribution functions corresponding to the last two 
variants will be constructed in the following subsection. The hard spheres 
model (H.S.) and Maxwell's molecules (M.M.) will be considered. 



5 Distribution functions of the second quasiequilibrium approxi- 
mation for scattering rates 



5.1 First five moments and collision stress tensor 

Elsewhere below the local equilibrium /^(34) is chosen as the first quasi- 
equilibrium approximation. 

Let us choose v ik = mViV k (36) as the microscopic density /i(v) of the second 
quasiequilibrium state. Let us write down the corresponding scattering rate 
(collision stress tensor) A ik in the form, 

A ifc = -(</?, v Qik ), (70) 
where 

VQik(v) = m J w(v',v[ | v 1 ,v)f (0) (v 1 ) 

x |(ui-«i)(ujfc - u k ) - ^5 ik (v - u,) 2 j dv'dv' 1 dv 1 (71) 
is the microscopic density of the scattering rate A ik . 

The quasiequilibrium distribution function of the second quasiequilibrium ap- 
proximation for fixed scattering rates (70) is determined as the solution to the 
problem 

(</?, ip) — > min for 

(<P,VQik) = -A ifc . (72) 
The method of Lagrange multipliers yields 
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(f(v) = \ikVQik(v), 
Kk(vQik, VQls) = A; s , 



(73) 



where are the Lagrange multipliers. 

In the examples of collision models considered below (and in general, for cen- 
trally symmetric interactions) vq^ is of the form 

u Qlk (v) = (i - P(°))i/ tt (»)$((« - uf), (74) 



where (1 — P^)u ik is determined by relationship (37) only and function $ 
depends only on the absolute value of the peculiar velocity (v — u). Then 

\k — f^ik'i 

r- 1 = (2/15) ($ 2 ((v - u) 2 ), (v - u) 4 ) , (75) 

and the distribution function of the second quasiequilibrium approximation 
for scattering rates (70) is given by the expression 

f = f(°\l + rA lkf i Qlk ). (76) 

The form of the function — it) 2 ), and the value of the parameter r are 
determined by the model of particle's interaction. In the Appendix, they are 
found for hard spheres and Maxwell molecules models (see (148)-(153)). The 
distribution function (76) is given by the following expressions: 

For Maxwell molecules: 



/ = / 



(o) 



X 



1 + ^ M m(2P 2 k B T)- 1 A ik {[ Vl - u t )(v k - u k ) - l -8 ik {v - u) 2 ^ } , 



m.m. k B TV2m 



to ' ' ~ 3^(5)^' (77) 



where fj^ M ' is viscosity coefficient in the first approximation of the Chapman- 
Enskog method (it is exact in the case of Maxwell molecules), k is a force 
constant, A 2 (5) is a number, A 2 (5) 0.436 (see [34]); 

For the hard spheres model: 
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/ = / (0) 




y 2 )(i + y 2 ) 



S ik (v 




^• s - = {5yJk B Tm)/{16y/^a 2 ), 



where r is a number represented as follows: 



^ = yJ /«- 11/2 /3(y)/3(^)7(y)7(^) 



+1 +1 



x (16a 2 + 28a(7(y) + 7(2)) + 63-f(y)j(z)) dy dz, 
a = l + y 2 + z 2 , (3(y) = l + y 2 , 1 (y) = \-y\ 



(79) 



Numerical value of r 1 is 5.212, to third decimal point accuracy. 

In the mixed description, the distribution function of the second quasiequilib- 
rium approximation under fixed values of the moments and of the scattering 
rates corresponding to the microscopic density (36) is determined as a solution 
of the problem 

(ip, ip) — > min for (80) 



((i - P {0) )v ik ,Lp) = a ik , 
("Qik, V?) = A ife . 

Taking into account the relation (74), we obtain the solution of the problem 
(80) in the form, 

(p(v) = (\ lk <5>((v - u) 2 ) + p ik )(( Vi - Ui)(v k - u k ) - (1/S)5 ik (v - u) 2 ).(81) 



Lagrange multipliers \ k ,(3i k are determined from the system of linear equa- 
tions, 

ms~ l \ ik + 2Pk B Tm~ 1 (3 ik = a ik , 



mr 1 \ ik + ms l (3 ik 



= A, 



(82) 
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where 

s- 1 = (2/15)($((w - u) 2 ), (v - u) 4 ). 



(83) 



If the solvability condition of the system (82) is satisfied, 

D = m 2 s- 2 - 2Pk B Tr~ 1 ^ 0, (84) 

then the distribution function of the second quasiequilibrium approximation 
exists and takes the form 

f = f(°) {l + ( m V 2 -2Pk B Tr- 1 )- 1 (85) 
x [(ms-^ik - 2Pk B Tm- 1 A tk )<S>((v - u) 2 ) 

+ (ms -1 A ifc - mr _1 (jj fe )]((^ - Ui)(v k - u k ) - (l/3)5 ik (v - u) 2 )} . 



The condition (84) means independence of the set of moments <Ji k from the 
scattering rates Ai k . If this condition is not satisfied, then the scattering rates 
Aik can be represented in the form of linear combinations of Oi k (with coeffi- 
cients depending on the hydrodynamic moments). Then the closed by means 
of (76) equations of the second chain are equivalent to the ten moment Grad 
equations, while the mixed chain does not exist. This happens only in the case 
of Maxwell molecules. Indeed, in this case 

s- 1 = 2P 2 k B T(m 2 ^ M -y 1 ;D = 0. 

The transformation changing A ik to a ik is 

/^•A^P- 1 = a ik . (86) 



For hard spheres: 

+i 



_, hP 2 k B T 
s - 



§-\ ~s-' = f i(y)(P(y)r 7/2 (f3(y) + li(y)) dy.(87) 



-1 



The numerical value of s 1 is 1.115 to third decimal point. The condition (83) 
takes the form, 



„ 25 P 2 k B T\ 2 . 2 , 
32 \mfif s - J y ' 



(88) 
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Consequently, for the hard spheres model the distribution function of the 
second quasiequilibrium approximation of the mixed chain exists and is de- 
termined by the expression 



/ = /(°) {l + m(APk B T(:s- 2 - r 1 ))- 1 



x 



-1 Sv^h 



5P 



t$- s -A ik \ J exp I - 



m{v-uf 2 N 
2k B T y , 



x(l-y 2 )(l+y 2 )( m( l y (l-y 2 )+2)dy 



+2 5' 



"I 8 ^ ..H.S. A 



2k B T 



x((vi - Ui)(v k - u k ) - ^5 ik (v - w ) 2 )} 



(89) 



5.2 First five moments, collision stress tensor, and collision heat flux vector 



Distribution function of the second quasiequilibrium approximation which 
takes into account the collision heat flux vector Q is constructed in a sim- 
ilar way. The microscopic density £gj is 

Coi(v) = J w(v',v[ | v,v 1 )f^(v 1 ){(i-P^)^ dv'dvldv!. (90) 



The desired distribution functions are the solutions to the following optimiza- 
tion problems: for the second chain it is the solution to the problem (72) with 
the additional constraints, 



m((p,£ Q i) = Qi- 



(91) 



For the mixed chain, the distribution functions is the solution to the problem 
(80) with additional conditions, 



m((p,^ Q i) = Qi, 
m(^(i-PW)&) = ft. 



(92) 
(93) 



Here £j = ViV 2 /2 (see (39)). In the Appendix functions are found for 
Maxwell molecules and hard spheres (see (153)- (158)). Since 
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(tQi,VQk j ) = ((i-P l ° ) )Zi,VQk j ) 

= fa, (i - p (0) K) = ((i - P {0 %, (i - P (0) > k3 ) = o, (94) 



the conditions (91) are linearly independent from the constraints of the prob- 
lem (72), and the conditions (93) do not depend on the constraints of the 
problem (80). 

Distribution function of the second quasiequilibrium approximation of the 
second chain for fixed A ik , Qi is of the form, 

/ = + rA tk u Qlk + riQitoi). (95) 



The parameter 77 is determined by the relation 

ir 1 = (l/3)(fo,fci)- (96) 

According to (157), for Maxwell molecules 

„ = ^ L (97) 

and the distribution function (95) is 
/ = / (0) 

x {l + l if M -m{2P 2 k B T)- 1 \ k {{v l - u t )(v k - u k ) - (l/3)5 lk (v - uf) 
+ ^MP^T)-^ - «0 - l) } • (98) 

For hard spheres (see Appendix) 
_ , 64m 3 (/i"- s -) 2 

V ~ ^125P 3 (A;bT) 2 ' ( " } 

where 77 is a number equal to 16.077 to third decimal point accuracy. 
The distribution function (95) for hard spheres takes the form 

{ 2\/2frau?' S- y / m(v — u) 2 \ „ s „ , 
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X 



' m(v — u) 2 
v 2k B T 
2v / 2r/m 3 /XQ- s 



7(y) + 2^J dy ((v* - Ui){v k - u k ) - ^5 ik (v - u) 2 ^j 



25P 2 (k B T) 2 
+1 
x J exp 

-l 



■Qi 



[Vi - Ui) \ {v-u) 



m(v — u)' 
2k B T 



y 2 PivMv) 



2 5k B T 
m 



2k B T 



■7(2/) + 2 dy 



- ^) (w - u) 2 J exp ^- 



2k B T 



(100) 



The functions /3(y),'y(y),cr(y) and are 

/%) = 1 + y 2 , 7 (y) = 1 - 2/ 2 , <r(j/) = 2/ 2 (l - 2/ 2 ), %) = 3y 2 - 1. (101) 



The condition of existence of the second quasiequilibrium approximation of 
the mixed chain (84) should be supplemented with the requirement 

R = m 2 T- 2 - hP ^ T )\ -i ± o. (102) 
2m 



Here 



(103) 



For Maxwell molecules 

r- 1 = (5P 2 k 2 B T 2 )/(3^ M -m 3 ), 

and the solvability condition (102) is not satisfied. Distribution function of 
the second quasiequilibrium approximation of mixed chain does not exist for 
Maxwell molecules. The variables Qi are changed to q± by the transformation 

3^ M -Q t = 2P gi . (104) 



For hard spheres, 



= T- 1 



2h{Pk B Tf 



(105) 



where 
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1 



-1 



r- 1 = ^ j /r 9/2 (y)7(y){63( 7 (y) + a(y)) 

-l 

+7(3(y)(4-10 1 (y) + 25(y)-5a(y)) 
+(3 2 (y)(25 1 (y) - 10%) - 40) + 20(3 3 (y)} dy. 



(106) 



The numerical value of r 1 is about 4.322. Then the condition (102) is verified: 

R « 66m~ 4 (PA;BT) 4 (/i^ s -) 2 . 

Finally, for the fixed values of cr^, A^, and <5« the distribution function 
of the second quasiequilibrium approximation of the second chain for hard 
spheres is of the form, 



/ = /(«> 1 + 



m 



4PkuT 



x 



+1 



5P 



m(v — u) 9 

-y 



2A; B T 



xP(yh(y) [ mi 2 kB T )2 l{y) + 2 ) d ^ + 2 ( 5 

x ({vi - Ui){v k - u fc ) - ^5 ik {v - u) 



"~ 1 ^^o" S Afc-r- 1 <T ifc 



m 



(f-^-r 1 ) 



10(Pfc B T) 2 

X [ («i - Mi) ((«-«) 



T 9* - 



5P 



5/cbP 
m 



)/exp(- 



A*0 



2/crT ^ , 



xP(vh(y) \, ' i{y) + 2 dy + - - u) s 



x 



+1 / 
/exp ^- 



v 2fc B T 



2k B T 



m{v — u) 
2k^T~ 



5A; R T S 



5P 



m 



(107) 



Thus, the expressions (77), (78), (89), (98), (100) and (107) give distribution 
functions of the second quasiequilibrium approximation of the second and 
mixed hydrodynamic chains for Maxwell molecules and hard spheres. They 
are analogues of ten- and thirteen-moment Grad approximations (38), (42). 

The next step is to close the second and mixed hydrodynamic chains by means 
of the found distribution functions. 
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6 Closure of the second and mixed hydro dynamic chains 



6.1 Second chain, Maxwell molecules 



The distribution function of the second quasiequilibrium approximation under 
fixed A ik for Maxwell molecules (77) presents the simplest example of the 
closure of the first (47) and second (69) hydrodynamic chains. With the help 
of it, we obtain from (47) the following transport equations for the moments 
of the first (local equilibrium) approximation: 

dtp + di(uip) = 0; 

p(d t u k + uAu k ) + d k P + ^(P-Vo^-A^) = 0; 

^(d t P + Ui diP) + ^Pdim + p- l pf M -A ik d iUk = 0. (108) 

Now, let us from the scattering rate transport chain (69) find an equation for 
A ik which closes the system (70). Substituting (77) into (69), we obtain after 
some computation: 



2 

d t A ik + d s (u s A ik ) + A is d s u k + A ks d s Ui - -8 ik A ls d s ui 

+p 2 (pf M -)- 1 (d iUk + d kUl - h tk d sUl 

+P(p^ M T 1 \k + A tk d s u s = 0. (109) 

For comparison, let us give ten-moment Grad equations obtained when closing 
the chain (46) by the distribution functions (38): 



dtp + di(uip) = 0; 

p{d t u k + UidiU k ) + d k P + diU ik = 0; 

-(d t P + mdiP) + ^Pd lUl + aikdiUk = 0; (110) 

dt^ik + d s (u s o- ik ) + P (diU k + d k ui - ^6 ik d s u}j 

2 

+a is d s u k + a ks d s Ui - -5 ik ai s d s ui + P(pf M -y 1 a ik = 0. (Ill) 

Using the explicit form of p% LM -(77), it is easy to verify that the transforma- 
tion (86) maps the systems (108), (109) and (110) into one another. This is a 
consequence of the degeneration of the mixed hydrodynamic chain which was 
already discussed. The systems (108), (109) and (110) are essentially equiva- 
lent. These specific properties of Maxwell molecules result from the fact that 
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for them the microscopic densities (1 — P^)viV k and (i — P^)viV 2 are eigen 
functions of the linearized collision integral. 

6.2 Second chain, hard spheres 



We now turn our attention to the closure of the second and of the mixed 
hydrodynamic chains for the hard spheres model. Substituting the distribution 
function (78) into (46) and (69), we obtain an analogue of the systems (108) 
and (109) (second chain, hard spheres): 



d t p + d i (u i p) = 0; (112) 
p(d t u k + uAu k ) + d k P + fsT 1 • -^dM^P- 1 ^) = 0; 



d t A ik + d s (u s A ik ) + fa 1 (d s u s )A ik + H .s. [diU k + d k Ui - -5 ik d s u s 

2 



/ 2 
+f(fii + a 2 ) f A is d s u k + A ks d s Ui - -5 ifc A is 9 s 

+f(ai + a 3 ) ^A is i9 fe M s + A ks diU s - -5 ik A ls d s u^ + (Pra / Ho' S ')A ik = 0. 

The dimensionless parameters do, &i, 5.2 and a 3 are determined by the quadra- 
tures 



+1+1 

Si = -/ / /%)/W(*)7(v)«- w/2 (v,*) 
-i -i 

x{99 7 (y) 7 (z)( 7 (z) - 1) + 18c*(y,z)(2 7 (z)( 7 (z) - 1) 

+4 7 (y)(4 7 (z) - 3)) + 8a 2 (y,z)(^(z) - 3)}dydz; 
+1+1 

~ a2 = YeI I ^W7(i/)7 2 W^ 11/2 (y^){63 7 (y) 7 (z) 

-1 -1 

+Ua(y, z)(3j(y) + 2 7 (z)) + 24a%, z)} dy dz; 
+1+1 

5 3 = ^/ / «- 11/2 (y,^(i/)/3(^) 7 (y) 7 (^) 
-i -i 

x^O/h^)^) - 1) + 14(2 7 (z)( 7 (z) - 1) 

+ 7 (y)(3 7 (z) - 2))a(j/, z) + 8a 2 (y, z)(3 7 (*) - 2)}dj/^; (113) 
+1+1+1 

a ° w 15367! / / I mx ' y ' 
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X7(x)7(y)7(^){10395 7 (a:)7(|/)7(^) + 3780^(x, y, z) 

x (7(^)7(2/) + l{x)l{z) + i{y)i{z)) + 1680^ 2 (x, y, *) 

x(7(x) +7(1/) +7(2;)) + 960^ 3 (a:,?/,2;)}da;d|/d2;; 

ip(x,y,z) = l + x 2 + y 2 + z 2 . (114) 

Their numerical values are a>i m 0.36, a 2 ~ 5.59, a 3 ps 0.38, a ~ 2.92 to 
second decimal point. 



6.3 Mixed chain 



The closure of the mixed hydrodynamic chain with the functions (89) gives 
the following modification of the system of equations (113): 



dtp + diiuip) = 0; 

p{d t u k + UidiU k ) + d k P + diO ik = 0; 
3 5 

-(<%P + UidiP) + -PdiUi + a ik diU k = 0; 
d t cr ik + d s (u s (j ik ) + P (diU k + <9 fc Wi - ^5 ik d s u s ^j 
+a is d s u k + a ks d s Ui - -5 ik ai s d s ui + A ik = 0; 
<9 t A ifc + d s (u s A ik ) + g 8v /2^ H .s. (^ iUfc + dkUi " ^ikdsUsj 

5P f«l/Q \ 

+ I75F(F^ I )i^ ( ''" s)a, ' = 

1 / 2 

+ -(a x + a 2 ) ( a is d s u k + a ks d s Ui - -5 ik ai s d s ui 

+±(a 1 + as) (a iS d k u s + a k Au s - 2 -5 ik a lsds u^ 

~ ( 2 
+b 1 (d s u s )p ik + (61 + 6 2 ) U3 is d s u k + l3 ks d s Ui - -5 ik f3i s d s ui 

~ ( 2 
+(bi + h) U3 is d k u s + p ks diU s - -5 ik (3i s d s ui 

+ 872 W ^)'(s- 2 - Fj I WS^' + a ° ait ) = 0: (115) 



8 \/2 H S a 

«ifc = s (Jjfe • /i ' "Ajfc; 

„ 1 8 \/2 H S a 1 



(116) 



It is clear from the analysis of distribution functions of the second quasiequi- 
librium approximations of the second hydrodynamic chain that in the Grad 
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moment method the function $(c 2 ) is substituted by a constant. Finally, let 
us note the simplest consequence of the variability of function $(c 2 ). If /j,q is 
multiplied with a small parameter (Knudsen number Kn equal to the ratio of 
the main free path to the characteristic spatial scale of variations of hydro- 
dynamic values), then the first with respect to Kn approximation of collision 
stress tensor A-^ has the form, 

A® = P (d iUk + d k u t - ^8 ik d s v}j (117) 



for Maxwell molecules, and 



(o) 5f 



A-*! = -^j^-P ( d iUk + d k u t - -5 tk 8 s u s ) (118) 



for hard spheres. Substitution of these expressions into the momentum equa- 
tions results in the Navier-Stokes equations with effective viscosity coefficients 

fJ>eB, 

V* = (119) 
for Maxwell molecules and 



for hard spheres. When using ten-moment Grad approximation which does 
not distinguish Maxwell molecules and hard spheres, we obtain /x e jr = f$' S - 
Some consequences of this fact are studied below in Sect. 7. 



7 Alternative Grad equations and a "new determination of molec- 
ular dimensions" (revisited) 



7.1 Nonlinear junctionals instead of moments in the closure problem 



Here we apply the method developed in the previous sections to a classical 
problem: determination of molecular dimensions (as diameters of equivalent 
hard spheres) from experimental viscosity data. Scattering rates (moments 
of collision integral) are treated as new independent variables, and as an al- 
ternative to moments of the distribution function, to describe the rarefied 
gas near local equilibrium. A version of entropy maximum principle is used 
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to derive the Grad-like description in terms of a finite number of scattering 
rates. New equations are compared to the Grad moment system in the heat 
non-conductive case. Estimations for hard spheres demonstrate, in particular, 
some 10% excess of the viscosity coefficient resulting from the scattering rate 
description, as compared to the Grad moment estimation. All necessary de- 
tails of the second chain formalism that are necessary for this example are 
explained below. 

Here we consider a new system of non-hydrodynamic variables, scattering rates 
M Q {f): 

M Qtll2l3 (f)= J n ili2i3 Q(f,f)dv; (121) 

which, by definition, are the moments of the Boltzmann collision integral 
Q(fJ): 



Q(fJ) = J Wi, W, Wl ) {/(t/)/(«i) -f(v)f( Vl )} dv'dvldVL 

Here w is the probability density of a change of the velocities, (v, V\) — > 
(V, v'i), of the two particles after their encounter, and w is defined by a model 
of pair interactions. The description in terms of the scattering rates Mq (121) 
is alternative to the usually treated description in terms of the moments M: 
M hi2hU) = f Vhhisfdv. 

A reason to consider scattering rates instead of the moments is that Mq 
(121) reflect features of the interactions because of the w incorporated in their 
definition, while the moments do not. For this reason we can expect that, in 
general, a description with a finite number of scattering rates will be more 
informative than a description provided by the same number of their moment 
counterparts. 

To come to the Grad-like equations in terms of the scattering rates, we have 
to complete the following two steps: (i) to derive a hierarchy of transport 
equations for p, u, P, and Mg ili2i3 in a neighborhood of the local Maxwell 
states fo(p,u,P); (ii) to truncate this hierarchy, and to come to a closed set 
of equations with respect to p, u, P, and a finite number of scattering rates. 

In the step (i), we derive a description with infinite number of variables, which 
is formally equivalent both to the Boltzmann equation near the local equilib- 
rium, and to the description with an infinite number of moments. The ap- 
proximation comes into play in the step (ii) where we reduce the description 
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to a finite number of variables. The difference between the moment and the 
alternative description occurs at this point. 

The program of two steps, (i) and (ii), is similar to what is done in the Grad 
method [20], with the only exception (and this is important) that we should 
always use scattering rates as independent variables and not to expand them 
into series in moments. Consequently, we use a method of a closure in the step 2 
that does not refer to the moment expansions. Major steps of the computation 
will be presented below. 

7.2 Linearization 

To complete the step (i), we represent / as /o(l + ip), where fo is the local 
Maxwellian, and we linearize the scattering rates (121) with respect to ip: 



Here Lq is the usual linearized collision integral, divided by / . Though AMq 
are linear in ip, they are not moments because their microscopic densities, 
A/iQ, are not velocity polynomials for a general case of w. 

It is not difficult to derive the corresponding hierarchy of transport equations 
for variables AMq ili2i3 , p, u, and P (we refer further to this hierarchy as to the 
alternative chain): one has to calculate the time derivative of the scattering 
rates (121) due to the Boltzmann equation, in the linear approximation (122), 
and to complete the system with the five known balance equations for the 
hydrodynamic moments (scattering rates of the hydrodynamic moments are 
equal to zero due to conservation laws). The structure of the alternative chain 
is quite similar to that of the usual moment transport chain, and for this reason 
we do not reproduce it here (details of calculations can be found in [45]). One 
should only keep in mind that the stress tensor and the heat flux vector in the 
balance equations for u and P are no more independent variables, and they 
are expressed in terms of AMg ili2i3 , p, u, and P. 

7.3 Truncating the chain 

To truncate the alternative chain (step (ii)), we have, first, to choose a finite set 
of "essential" scattering rates (122), and, second, to obtain the distribution 
functions which depend parametrically only on p, u, P, and on the chosen 




(122) 




(^1112*3) ■ 
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set of scattering rates. We will restrict our consideration to a single non- 
hydrodynamic variable, CQij, which is the counterpart of the stress tensor a^. 
This choice corresponds to the polynomial mViVj in the expressions (121) and 
(122), and the resulting equations will be alternative to the 10 moment Grad 
system 2 . For a spherically symmetric interaction, the expression for gq^ may 
be written: 



(T Qij( ( f) = J An Qij f (pdv; (123) 

P r 1 i 

Afi Qij = LQimViVj) = j^S Q (c 2 ) \ c iCj - -5 ijC 2 \ . 

Vqo{J- ) L 3 J 



Here 7]qq(T) is the first Sonine polynomial approximation of the Chapman- 
Enskog viscosity coefficient (VC) [34], and, as usual, c = \J^§f{v — it). The 
scalar dimensionless function Sq depends only on c 2 , and its form depends on 
the choice of interaction w. 



7.4 Entropy maximization 



Next, we find the functions 

f*(p,u,P,a Qij ) = f (p,u,P)(l + ip*(p,u,P,a Qij )) 

which maximize the Boltzmann entropy S(f) in a neighborhood of fo (the 
quadratic approximation to the entropy is valid within the accuracy of our 
consideration), for fixed values of ctq^-. That is, ip* is a solution to the following 
conditional variational problem: 

AS(p) = - y / /oy? 2 dv - max, (124) 
i) J A/j, Qij f pdv = a Qij ; ii) J {l, v, w 2 } f p dv = 0. 

The second (homogeneous) condition in (124) reflects that a deviation p from 
the state f is due only to non-hydrodynamic degrees of freedom, and it is 
straightforwardly satisfied for A/iQij (123). 

Notice, that if we turn to the usual moment description, then condition (i) 
in (124) would fix the stress tensor instead of its scattering counterpart 

2 To get the alternative to the 13 moment Grad equations, one should take into 
account the scattering counterpart of the heat flux, qqi = m J Vi^Q(f, f) dv. 
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UQij. Then the resulting function f*(p,u, P,&ij) will be exactly the 10 mo- 
ment Grad approximation. It can be shown that a choice of any finite set 
of higher moments as the constraint (i) in (124) results in the corresponding 
Grad approximation. In that sense our method of constructing /* is a direct 
generalization of the Grad method onto the alternative description. 

The Lagrange multipliers method gives straightforwardly the solution to the 
problem (124). After the alternative chain is closed with the functions f*(p, u, P, (TQij), 
the step (ii) is completed, and we arrive at a set of equations with respect to 
the variables p, u, P, and ctq^-. Switching to the variables Qj = n^aQij, we 
have: 



dtn + di(nui) = 0; (125) 
p(d t u k + uAu k ) + d k P + d t j ^p^ Cifcj = 0; (126) 

^(d t P + Ul d t P) + \pd iUi + j ^p% | d iUk = 0; (127) 

2 

dtCik + d s (u s ( ik ) + {CksdsUi + ( is d s u k - -5i k Crsd s u r } (128) 

+ 1 1Q - — - 1 dkdsUs ^7=—(diU k + d k Ui - ^S ik d s u s ) 

{ r Q J VQo{T)n 3 

rCifc — U. 



rqVQo(T) 



Here dt = d/dt,di = d/dxi, summation in two repeated indices is assumed, 
and the coefficients tq, [3q, and «q are defined with the help of the function 
Sq (123) as follows: 



v 

oo 

olq = ^ / e^c e S Q (c 2 )R Q (c 2 ) dc. (129) 

* 

The function Rq(c 2 ) in the last expression is defined due to the action of the 
operator Lq on the function Sq(c 2 )(ciCj — |5yC 2 ): 

PI 1 

RQ(c 2 )(ciCj - -Sijc 2 ) = L Q (S Q (c 2 )(c i c j - -SijC 2 )). (130) 

Vqo <j <j 
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Finally, the parameter 7g in (125-129) reflects the temperature dependence of 
the VC: 

lQ A Vqo(T)[ dT ))■ 

The set of ten equations (125-129) is alternative to the 10 moment Grad 
equations. 

7.5 A new determination of molecular dimensions (revisited) 

The observation already made is that for Maxwell molecules we have: 5* M ' M ' = 

1, and r^ M - oc T; thus 7 M - M - = (3 MM - = 0, r MM - = a MM - = \, and (125-129) 

becomes the 10 moment Grad system under a simple change of variables XQj = 

o~ij, where A is the proportionality coefficient in the temperature dependence 
oir] M.M. 

These properties (the function Sq is a constant, and the VC is proportional to 
T) are true only for Maxwell molecules. For all other interactions, the function 
Sq is not identical to one, and the VC t]q (T) is not proportional to T. Thus, 
the shortened alternative description is not equivalent indeed to the Grad 
moment description. In particular, for hard spheres, the exact expression for 
the function S n ' S ' (123) reads: 

S R - S - = ^ / exp(-c¥)(l - t 4 ) (c 2 (l - t 2 ) + 2) df; (131) 
o 

% H - S - oc VT. 

Thus, 7 H ' S ' = |, and ~ 0.07, and the equation for the function (129) 
contains a nonlinear term, 

9 U - S -Q k d s u s , (132) 

where 6 |H S - w 0.19. This term is missing in the Grad 10 moment equation. 

Finally, let us evaluate the VC which results from the alternative description 
(125-129). Following Grad's arguments [20], we see that, if the relaxation of 
( ik is fast compared to the hydrodynamic variables, then the two last terms in 
the equation for Qk (125-129) become dominant, and the equation for u casts 
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Fig. 3. Approximations for hard spheres: bold line - function S* H ' S ', solid line 
approximation S^' S ', dotted line - Grad moment approximation. 

into the standard Navier-Stokes form with an effective VC i]q cS : 



1 



2«q 



Vqo- 



(133) 



For Maxwell molecules, we easily derive that the coefficient ckq in (133) is equal 
to |. Thus, as one expects, the effective VC (133) is equal to the Grad value, 
which, in turn, is equal to the exact value in the frames of the Chapman- 
Enskog method for this model. 

For all interactions different from the Maxwell molecules, the VC i]q c s (133) 
is not equal to r]Q . For hard spheres, in particular, a computation of the VC 
(133) requires information about the function R us - (130). This is achieved 
upon a substitution of the function S R ' S ' (131) into (130). Further, we have to 
compute the action of the operator L H S ' on the function S (ciCj — |<5jjC 2 ), 
which is rather complicated. However, the VC can be relatively easily 
estimated by using a function S^ s - = ^75 (1 + 7C 2 ), instead of the function 
gH.s.^ m (130). Indeed, the function S^" S ' is tangent to the function S K ' S- 
at c 2 = 0, and is its majorant (see Fig. 3). Substituting S^ s - into (130), 
and computing the action of the collision integral, we find the approximation 
R^' S ", thereafter we evaluate the integral a H ' s ' (129), and finally come to the 
following expression: 



n H.s. > 75264 h S 
VeS ~ 67237^ 



1.12 % H - S - 



(134) 



Thus, for hard spheres, the description in terms of scattering rates results in 
the VC of more than 10% higher than in the Grad moment description. 



A discussion of the results concerns the following two items. 
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Table 1 

Three virial coefficients: experimental B exp , classical Bq [46], and reduced B e g for 
three gases at T = 500K 







Bo 


B c ft 


Argon 


8.4 


60.9 


50.5 


Helium 


10.8 


21.9 


18.2 


Nitrogen 


168 


66.5 


55.2 



1. Having two not equivalent descriptions which were obtained within one 
method, we may ask: which is more relevant? A simple test is to compare char- 
acteristic times of an approach to hydrodynamic regime. We have tq ~ Vo' s '/P 
for 10-moment description, and r a ~ i]^ s -/P for alternative description. As 
r a > tq, we see that scattering rate decay slower than corresponding moment, 
hence, at least for rigid spheres, the alternative description is more relevant. 
For Maxwell molecules both the descriptions are, of course, equivalent. 

2. The VC r]fg S - (134) has the same temperature dependence as r^' 3 ', and also 
the same dependence on a scaling parameter (a diameter of the sphere). In the 
classical book [34] (pp. 228-229), "sizes" of molecules are presented, assuming 
that a molecule is represented with an equivalent sphere and VC is estimated 
as 7]^- s -. Since our estimation of VC differs only by a dimensionless factor from 
7/J 1 ' 3 ', it is straightforward to conclude that effective sizes of molecules will be 
reduced by the factor b, where 

b = y/ri-s-frfe 8 - » 0.94. 

Further, it is well known that sizes of molecules estimated via viscosity in [34] 
disagree with the estimation via the virial expansion of the equation of state. 
In particular, in book [46], p. 5, the measured second virial coefficient B exp 
was compared with the calculated B , in which the diameter of the sphere was 
taken from the viscosity data. The reduction of the diameter by factor b gives 
B c fi = b 3 B . The values B exp and B [46] are compared with B e s in the Table 
1 for three gases at T = 500K. The results for argon and helium are better 
for B e s, while for nitrogen B e s is worth than B . However, both B and B c s 
are far from the experimental values. 

Hard spheres is, of course, an oversimplified model of interaction, and the 
comparison presented does not allow for a decision between rj^ 3 ' and r^ 3 '. 
However, this simple example illustrates to what extend the correction to the 
VC can affect a comparison with experiment. Indeed, as it is well known, the 
first-order Sonine polynomial computation for the Lennard- Jones (LJ) poten- 
tial gives a very good fit of the temperature dependence of the VC for all 
noble gases [47], subject to a proper choice of the two unknown scaling pa- 
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rameters of the LJ potential 3 . We may expect that a dimensionless correction 
of the VC for the LJ potential might be of the same order as above for rigid 
spheres. However, the functional character of the temperature dependence will 
not be affected, and a fit will be obtained subject to a different choice of the 
molecular parameters of the LJ potential. The five-parametric family of pair 
potentials was discussed in Ref . [48] . These five constants for each pair poten- 
tial have been determined by a fit to experimental data with some additional 
input from theory. After that, the Chapman-Enskog formulas for the second 
virial coefficient and main transport coefficients give satisfactory description 
of experimental data [48]. Such a semiempirical approach that combines fitting 
with kinetic theory might be very successful in experimental data description, 
but does not allow us to make a choice between hierarchies. This choice re- 
quires less flexibility in the potential construction. The best solution here is 
independent determination of the interaction potential without references to 
transport coefficients or thermodynamic data. 



Appendix: Formulas of the second quasiequilibrium approximation 
of the second and mixed hydrodynamic chains for Maxwell molecules 
and hard spheres 



Let us write uq ik (71) in the standard form: 



VQik (135) 
= j / (0) I «i - v | - Ui)(v k - u k ) - -S ik (v - uf^bdbdedv!, 

where b is the impact parameter, e is the angle between the plane containing 
the trajectory of the particle being scattered in the system of the center of 
mass and the plane containing the entering asymptote, the trajectory, and a 
certain fixed direction. It is convenient to switch to the dimensionless velocity 
c: 

771 \ 

(vi-ui) (136) 



2knT 



and to the dimensionless relative velocity g: 




(137) 



3 A comparison of molecular parameters of the LJ potential, as derived from the 
viscosity data, to those obtained from independent sources, can be found elsewhere, 
e.g. in [34], p. 237. 
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After standard transformations and integration with respect to e (see [34]) we 
obtain in (136) 



vc*k = — vr~ 1/2 (138) 
m 

x j exp(-c^)v?S 2) (5() {[c u - Ci)(c lk - c k ) - -<5 ifc (ci - c) 2 ^j dci. 

Here 

V^S 2) = |(l-cos 2 x)|«i-«|&(x)|^|^, (139) 



and x is an angle between the vectors g and g' . 

(2) 

The dependence of ip\ on the vector g is determined by the choice of the 
model of particle's interaction. 

For Maxwell molecules, 

^>=(t)%(5), (140) 



where k is a force constant, A 2 (5) is a number, A 2 (5) ~ 0.436. 
For the model of hard spheres 

(2) V2aUk B T\ l/2 . 

<Pi =—5— ci-c , 141 

3 \ m J 

where a is diameter of the sphere modelling the particle. 

Substituting (140) and (141) into (139), we transform the latter to the form: 

for Maxwell molecules 



3P/2k\ 1 / 2 .., , 2 , ( d d 1. d 9\ MM/2N 

= 4^ (™) A2(5) exp( " c } " 3^^J T (c ); 

T M.M. (c2) = I exp (_ a .2 _ 2xk( , k) dx . (M2) 



for hard spheres 
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Pa 2 /A; B T\ 1/2 / d d I d d \ h. s 2 

v Q ik = -7= exp -c — q^ifco— 7T~ T c ); 

2y/2m\7rmJ \dcidc k 3 dc s dc s J 

T H - s -(c 2 ) = J \ x\ exp(-x 2 - 2x k c k ) dx. (143) 



It is an easy matter to perform integration in (142), the integral is equal to 

7T 3 /V 2 . 



Therefore for Maxwell molecules, 



3 /2k\ ' ( 1 \ 

v Qik = 2 n7r ( — J ^2(5) f («< - Ui)(v k - u k ) - -S ik (v - u) 2 \ . (144) 

The integral T H - S - in (143) can be transformed as follows: 

+1 

T H - s -(c 2 ) = 2tt + 7T I exp(c 2 (l - y 2 ))c 2 (l + y 2 ) dy. (145) 



Then for the model of hard spheres, 



r— . 1 k B T\ 3 ^ 2 ( 1 , 

v Qik = V2ima [ — — I I qc fc - -d ifc c 



+1 

x | exp(-cV)(l+y 2 )(l-y 2 )(c 2 (l-y 2 )+2)dy. (146) 

-1 



Let us note a useful relationship: 



+1 



d n T H - s -/d(c 2 ) n = 71 J exp(c 2 (l - y 2 )) 
-1 

x(l + y 2 )(l - y 2 ) n -\c 2 {\ - y 2 ) + n) dy, n > 1. (147) 



Use the expressions for the viscosity coefficient yU which are obtained in the 
first approximation of the Chapman-Enskog method: 

for Maxwell molecules, 

_/2mxV 2 k B T 

"° '[—J 3^(5)' (148) 
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for hard spheres, 

h.s. _ 5(fcBTm)^ 



Transformation of (144), (146) to the form of (74) gives the following functions 
<H(v uf): 

for Maxwell molecules, 

$ = P/zi^-M.. ( 150 ) 

for hard spheres 



x (1 + y 2 )(l - V 2 ) { m{ l~J ] \ ^ ~ V 2 ) + 2) dy. (151) 



The parameter r from (75) is: 
for Maxwell molecules: 

r=(m^ M -) 2 /(2P 3 k B Ty, (152) 

for hard spheres: 

64 fm/i^Y 

r = f V ,° >-. (153) 

25P 3 k B T K J 

The dimensionless parameter f is represented as follows: 



+1+1 

f" 1 = ^ / / oT 11 /^)/^)^)^) 
-1 -1 

x (16a 2 + 28a(7(y) + 7(2)) + 637(^)7(2)) dy ds. (154) 



Here and below the following notations are used: 

(3{y) = l + y\ 7 (y) = l-</ 2 , a = l + y 2 + z 2 . (155) 
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Numerical value of f~ l is 5.212 to third decimal point. 
The parameter (83) is: 
for Maxwell molecules 

s- 1 = (2P 2 k B T)/ (m 2 ^- M -) ; (156) 

for hard spheres 



8m 2 /ig ' 



s _1 = g" 1 L 2 .,h.s. ■ (157) 



The dimensionless parameter s 1 is of the form 

r 1 = / i(y)P' 7/2 (y) (p(y) + \i{y)) dy. (158) 

Numerical value of is 1.115 to third decimal point. 
The scattering rate density (90) is of the form, 

Z Qi = y /2^y Jf°\v 1 )\v 1 -v\l [ c i (c 2 -^bdbdedv 1 . (159) 



Standard transformation of the expression {cj(c 2 — 5/2)} and integration with 
respect to e change (159) to the form, 



fa (160) 
J exp(-c 2 )v9S 2) (3(c 2 -c 2 )(c li -c i )-(c 1 -c) 2 (c li + Q))dc 1 . 



2itm 

Further, using the expressions (140) and (141) for ipf \ we obtain: 
for Maxwell molecules: 

1 /2 

eOi = 4f— ) ^(5)exp(-c 2 )AT M - M -(c 2 ); (161) 



47 



for hard spheres: 



^ = ^£^exp(-c 2 )AT H - s -(c 2 ). (162) 



The operator Di is of the form 

\d_d__d_ 3 c d_d_ _ \^ d_d_ 

4 dci dc s dc s 2 s dc s dci 2 1 dc s dc s 



The operator acts on the function ?/>(c 2 ) as follows: 

*fr(*A)+cj(*m. (164) 



rf(c 2 ) 2 1 V 2/ 4 V rf ( c2 ) 2 d{c 2 f)' 

From (161), (162) we obtain: 
for Maxwell molecules: 

^o* = 3^sr( u < - u <) (( v - «) 2 - ^^") 5 ( 165 ) 

for hard spheres: 



+ 1 / / 

y exp {-—k 

+(«» - Ui)(v - uf 



+ 1 / / \2 \ / / x2 



y exp \ 2A 

-1 x 



w) 2 ' ~ ' - ' , .m[v — u) 



x / ex P ( — 2kBT y J P(vh(y) 2A;bT + %) I 



The functions cr(y), 5(y) are of the form 

*(y) = y 2 (i-y 2 ), 6( y ) = tf-i. (167) 

The parameter 77 from (96) is: 
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for Maxwell molecules: 



9m 3 (fi¥ M -) 2 

"= iopWr (168) 



for hard spheres: 

^64m 3 (/io' s -) 2 

77 = *> 125P3(fc B T)2 • (169) 
The dimensionless parameter 77 is of the form 



+1+1 

V-* = J j ^mi{y)l{^- U/2 ^{l{y)l{z) + a{y)a{z) 
-1 -1 

63 

+ cr(y)7(z) + a(z)^(y)) + — a(2 7 (2/) + 27(2) - 57(^)7(2;) 

+ 2(a(y) + <t(*)) + 7 (z)5(y) + 7(2/)^) + <Ks/)*(*) + ^)5{y)) 

7 25 
+ -a 2 (4 - 10 7 (y) - 107(2;)) + -jrtyMz) + 2<%) (170) 

+ 25(z) - 5a(y) - 5a(z) - \{i{z)5{y) + >y(y)5(z) + 6(y)5(z)) 

+ I a 3 (_ 20 + ^( 7 (y) + 7 (*)) - 5(%) + 5(z))) + ^a 4 } dy dz. 



Numerical value of fj 1 is 0.622 to second decimal point. 
Finally, from (165), (167) we obtain r _1 (103): 
for Maxwell molecules 

1 _ 5(Pk B T)\ 

3^- M m 3 ' 1 ' 

for hard spheres 



T =T 



,3„H.S 

+1 

-■1 I 

r 



1 = ^/r 9/2 (y)7(y){63( 7 (y) + ^)) 

-1 

+ 7/%) (4 - 107(2/) + 25(2/) - 5(7(2/)) + 20/3%) 
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+ (3 2 {y){2h 1 {y) - 10%) - 40)} dy » 4.322. 



(172) 



8 Conclusion and outlook 

We developed the Triangle Entropy Method (TEM) for model reduction and 
demonstrated how it works for the Boltzmann equation. Moments of the Boltz- 
mann collision integral, or scattering rates are treated as independent variables 
rather than as infinite moment series. Three classes of reduced models are con- 
structed. The models of the first class involve only moments of distribution 
functions, and coincide with those of the Grad method in the Maximum En- 
tropy version. The models of the second type involves only scattering rates. 
Finally, the mixed description models involve both the moments and the scat- 
tering rates. TEM allows us to obtain all the closure formulas in explicit form, 
not only for the Maxwell molecules (as it is usual), but for hard spheres also. 
We found the new Boltzmann-kinetics estimations for the equivalent hard 
spheres radius for gases. 

The main benefits from TEM are: 

(1) It constructs the closure as a solution of linear equations, and, therefore, 
often gives it in an explicit form; 

(2) It provides the thermodynamic properties of reduced models, at least, 
locally; 

(3) It admits nonlinear functionals as macroscopic variables, this possibility is 
important for creation of non-equilibrium thermodynamics of non-linear 
fluxes, reaction rates, scattering rates, etc. 

The following fields for future TEM applications are important: 

• Modelling of nonequilibrium processes in gases (Boltzmann kinetics and its 
generalisations); 

• Chemical kinetics models with reaction rates as independent variables; 

• Kinetics of complex media (non-Newtonian liquids, polymers, etc.) with the 
Fokker-Planck equation as the basic kinetic description. 

Of course, TEM is a new computational version of the entropy maximum 
method (MaxEnt) that is in wide use in the Extended irreversible thermody- 
namics (EIT) [44] after works of Jaynes [3], Kogan and Rosonoer [5,6,7] and 
Lewis [8]. 

MaxEnt methods obtain the second wind because the high interest to non- 
classical nonextensive entropies [49,50]. In this sense, the Fokker-Plank equa- 
tion seems to be a very attractive example for MaxEnt method application, 



50 



and, in particular, for the application of TEM. This classical equation admits 
a broad class of Lyapunov functions, including nonextensive entropies. 



The Fokker-Planck equation (FPE) is a familiar model in various problems of 
nonequilibrium statistical physics We consider the FPE of the form 



dW(x,t) d_< 
dt dx I 



d d 
WJ-U+—W 
ox ox 



(173) 



Here, W(x, t) is the probability density over the configuration space x at time 
t, while U(x) and D(x) are the potential and the positively semi-definite 
((y,Dy) > 0) diffusion matrix. 

It is known that for the Boltzmann equation there exists only one universal 
Lyapunov functional: the entropy (we do not distinguish functionals which 
are related to each other by monotonic transformations). For the FPE there 
exists a whole family of universal Lyapunov functionals. Let h(a) be a convex 
function of one variable a > 0, h"(a) > 0, 

Sh[W] = - j W cq (x)h 



W(x,t) 
W eq (x) 



dx. 



(174) 



The density of production of the generalized entropy Sh, cr h , is nonnegative 

~W(x,t) 



a h (x) = W eq (x)h" 



dW (*>%D±?£m>0. (175) 



W cq (x) \\dx W cq (x) dx W eq (x) 



The most important variants for the choice of h are: 

• h(a) = a In a, and Sh is the Boltzmann-Gibbs-Shannon entropy (in the 
Kullback form [51,52]), 

• h(a) = a In a — elna, e > 0, and is the maximal family of additive 
entropies [53,54,55] (these entropies are additive for the composition of in- 
dependent subsystems). 

• h(a) = ^fzf, and is the family of Tsallis entropies [49,50]. These entropies 
are not additive, but become additive after a nonlinear monotonous trans- 
formation. This property can serve as a definition of the Tsallis entropies in 
the class of generalized entropies (174) [55]. 

The MaxEnt closure approximations for the Fokker-Planck equation are de- 
veloped for the classical Boltzmann-Gibbs-Shannon entropy in the Kullback 
form and standard linear functionals (moments) as macroscopic variables 
[52,27,28]. Now we are in position to develop the whole family of approxima- 
tion in explicit form (due to TEM), for classical and non-classical entropies, 
for linear and non-linear macroscopic variables. 
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The crucial question is: where to stop? Is it possible to decide, is this partic- 
ular model from the hierarchy sufficiently accurate, or we need to go ahead? 
Without criteria for such a decision we have infinite number of theories. 

The residual estimates are possible: we can estimate the defect of invariance 
(see Fig. 1). If it is too big (in comparison with the full right hand side J), 
then we should switch to the next system of hierarchy. If it vanishes, we could 
try the previous system. Normally, it is impossible to find one reduced model 
for all regimes, but it is possible to change the model during simulation. 

There exists one more benefit from the hierarchy. For each model we have a 
correspondent approximate slow invariant manifold Q{, and a vector field of 
reduced dynamics Jj which is defined at points from and is tangent to this 
manifold. This structure gives a possibility to estimate not the whole defect 
of invariance, but a "partial defect" Aj = J i+1 — Ji. Usually it is sufficient to 
estimate this partial defect of invariance, that is, to check whether the current 
model is the approximate slow invariant manifold for the next model up to 
desired accuracy. Examples of these estimates and applications are presented 
in Refs. [16,17,30,39]. We propose to use the flexible technology of modeling 
with adaptive choice of the model from hierarchy. This approach could be 
discussed as intermediate one between the classical one-model calculations 
and the equation-free approach [57], for example. 

We construct the quasiequilibrium hierarchy of models for a system with en- 
tropy growth. These systems relax to equilibrium points. But most interesting 
application is modeling of open systems. It is possible to use obtained hier- 
archy of models for open systems just by adding flows under the assumption 
that the fast motion and slow manifold do not change due to the system open- 
ing. For example, we usually use the Navier-Stokes equation for systems with 
external flows that do not relax to equilibrium. If the external flows are fast 
and the perturbation of slow manifold is significant, then the correspondent 
perturbation theory [16,17] modifies the model for open system. 

The inertial manifold [14,58,59] is the manifold where the limit behaviour of 
the system is located; it exponentially attracts motions when t — > oo. For a 
closed system the equilibrium (one point) is already the inertial manifold. In 
the theory of inertial manifolds the estimates of inertial manifolds dimension 
for several classes of (open) systems are created and finiteness of this dimension 
is proved. Inertial manifolds could be considered as the lowest level of any hier- 
archy of slow manifolds. They belong to all the slow invariant manifolds of the 
hierarchy. In our construction we build the hierarchy of infinite-dimensional 
approximate slow manifolds for the Boltzmann equation and do not try to 
find the smallest invariant manifolds for open systems. 

And, finally, we should ask the question: what chain is better, could we prove 



52 



that the second hierarchy with scattering rates instead of usual moments is 
better than the standard Grad hierarchy? We cannot prove this exactly, but 
can only argue plausibly that the second hierarchy should lead to dynamic in- 
variance faster, than the first one, and support this point of view by examples. 
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